<!-- XHTML 1.0 Strict -->
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en">
<meta name="author" content="Rachid Touzani" />
<meta name="keywords" content=" css, dropdowns, dropdown menu, drop-down, menu, navigation, nav, horizontal, vertical left-to-right, vertical right-to-left, horizontal linear, horizontal upwards, cross browser, internet explorer, ie, firefox, safari, opera, browser, lwis" />
<meta name="description" content="Clean, standards-friendly, modular framework for dropdown menus" />
<link href="../css/dropdown.vertical.css" media="screen" rel="stylesheet" type="text/css" />
<link href="../css/default.ultimate.css" media="screen" rel="stylesheet" type="text/css" />
<head>
<!-- Title -->
<title>OFELI - Documentation</title>
<!-- Stylesheet -->
<link rel="stylesheet" type="text/css" href="../doxygen.css" />
<link rel="stylesheet" type="text/css" href="../tabs.css" />
</head>
<body>
<p align="center">
  <img src="../im/ofeli.gif" alt="OFELI's Logo" width="300" height="90">
</p>
<p align="center"><font color="#203F5B" size="4" face="verdana,Helvetica">
<b>An Object Oriented Finite Element Library</b></font></p>
<!-- Complete Page-->
<div id="main">
<div id="mainmenu">
 <div class="text">
<!-- Menu -->
 <div class="tabs">
  <ul class="tablist">
   <li><a href="../index.html"><span>Home</span></a>
    <li><a href="../overview.html"><span>Overview</span></a>
    <li><a href="../extlib.html"><span>External Libraries</span></a>
    <li><a href="../citing.html"><span>Citing OFELI</span></a>
    <li><a href="../fformats.html"><span>File Formats</span></a>
    <li><a href="../demos.html"><span>Demo Codes</span></a>
  </ul>
 </div>
<!-- Generated by Doxygen 1.8.11 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="index.html"><span>Main&#160;Page</span></a></li>
      <li><a href="modules.html"><span>Modules</span></a></li>
      <li><a href="namespaces.html"><span>Namespaces</span></a></li>
      <li class="current"><a href="annotated.html"><span>Classes</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
        <span class="left">
          <img id="MSearchSelect" src="search/mag_sel.png"
               onmouseover="return searchBox.OnSearchSelectShow()"
               onmouseout="return searchBox.OnSearchSelectHide()"
               alt=""/>
          <input type="text" id="MSearchField" value="Search" accesskey="S"
               onfocus="searchBox.OnSearchFieldFocus(true)" 
               onblur="searchBox.OnSearchFieldFocus(false)" 
               onkeyup="searchBox.OnSearchFieldChange(event)"/>
          </span><span class="right">
            <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
          </span>
        </div>
      </li>
    </ul>
  </div>
  <div id="navrow2" class="tabs2">
    <ul class="tablist">
      <li><a href="annotated.html"><span>Class&#160;List</span></a></li>
      <li><a href="classes.html"><span>Class&#160;Index</span></a></li>
      <li><a href="hierarchy.html"><span>Class&#160;Hierarchy</span></a></li>
      <li><a href="functions.html"><span>Class&#160;Members</span></a></li>
    </ul>
  </div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
      <div id="nav-sync" class="sync"></div>
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
$(document).ready(function(){initNavTree('classOFELI_1_1PETScMatrix.html','');});
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
     onmouseover="return searchBox.OnSearchSelectShow()"
     onmouseout="return searchBox.OnSearchSelectHide()"
     onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>

<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0" 
        name="MSearchResults" id="MSearchResults">
</iframe>
</div>

<div class="header">
  <div class="summary">
<a href="#pub-methods">Public Member Functions</a> &#124;
<a href="classOFELI_1_1PETScMatrix-members.html">List of all members</a>  </div>
  <div class="headertitle">
<div class="title">PETScMatrix&lt; T_ &gt; Class Template Reference<div class="ingroups"><a class="el" href="group__OFELI.html">OFELI</a><a class="el" href="group__OFELI.html">OFELI</a> &raquo;  &#124; <a class="el" href="group__VectMat.html">Vector and Matrix</a></div></div>  </div>
</div><!--header-->
<div class="contents">

<p>To handle matrices in sparse storage format using the Petsc library.  
 <a href="classOFELI_1_1PETScMatrix.html#details">More...</a></p>
<table class="memberdecls">
<tr class="heading"><td colspan="2"><h2 class="groupheader"><a name="pub-methods"></a>
Public Member Functions</h2></td></tr>
<tr class="memitem:a9708182d1255298f5eb71c68d7270625"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a9708182d1255298f5eb71c68d7270625">PETScMatrix</a> ()</td></tr>
<tr class="memdesc:a9708182d1255298f5eb71c68d7270625"><td class="mdescLeft">&#160;</td><td class="mdescRight">Default constructor.  <a href="#a9708182d1255298f5eb71c68d7270625">More...</a><br /></td></tr>
<tr class="separator:a9708182d1255298f5eb71c68d7270625"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ac4143fbd170eb4ee9ddb38c0c0b53614"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ac4143fbd170eb4ee9ddb38c0c0b53614">PETScMatrix</a> (size_t nr, size_t nc)</td></tr>
<tr class="memdesc:ac4143fbd170eb4ee9ddb38c0c0b53614"><td class="mdescLeft">&#160;</td><td class="mdescRight">Constructor that initializes current instance as a dense matrix.  <a href="#ac4143fbd170eb4ee9ddb38c0c0b53614">More...</a><br /></td></tr>
<tr class="separator:ac4143fbd170eb4ee9ddb38c0c0b53614"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:af521d98af518239d3362748df7d02ab6"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#af521d98af518239d3362748df7d02ab6">PETScMatrix</a> (size_t size)</td></tr>
<tr class="memdesc:af521d98af518239d3362748df7d02ab6"><td class="mdescLeft">&#160;</td><td class="mdescRight">Constructor that initializes current instance as a dense matrix.  <a href="#af521d98af518239d3362748df7d02ab6">More...</a><br /></td></tr>
<tr class="separator:af521d98af518239d3362748df7d02ab6"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a91e7a1b127844bd9b2846db0c928d4a9"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a91e7a1b127844bd9b2846db0c928d4a9">PETScMatrix</a> (<a class="el" href="classOFELI_1_1Mesh.html">Mesh</a> &amp;mesh, size_t dof=0)</td></tr>
<tr class="memdesc:a91e7a1b127844bd9b2846db0c928d4a9"><td class="mdescLeft">&#160;</td><td class="mdescRight">Constructor using a <a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance.  <a href="#a91e7a1b127844bd9b2846db0c928d4a9">More...</a><br /></td></tr>
<tr class="separator:a91e7a1b127844bd9b2846db0c928d4a9"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a5ec47afc10afd3a677382bbbe4f06819"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a5ec47afc10afd3a677382bbbe4f06819">PETScMatrix</a> (const vector&lt; std::pair&lt; size_t, size_t &gt; &gt; &amp;I, int opt=1)</td></tr>
<tr class="memdesc:a5ec47afc10afd3a677382bbbe4f06819"><td class="mdescLeft">&#160;</td><td class="mdescRight">Constructor for a square matrix using non zero row and column indices.  <a href="#a5ec47afc10afd3a677382bbbe4f06819">More...</a><br /></td></tr>
<tr class="separator:a5ec47afc10afd3a677382bbbe4f06819"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a57718720311c145909fdf29cdc13eabf"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a57718720311c145909fdf29cdc13eabf"></a>
&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a57718720311c145909fdf29cdc13eabf">PETScMatrix</a> (const <a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> &amp;m)</td></tr>
<tr class="memdesc:a57718720311c145909fdf29cdc13eabf"><td class="mdescLeft">&#160;</td><td class="mdescRight">Copy constructor. <br /></td></tr>
<tr class="separator:a57718720311c145909fdf29cdc13eabf"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a838b1147e06cd9f9eb3795437d331471"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a838b1147e06cd9f9eb3795437d331471"></a>
&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a838b1147e06cd9f9eb3795437d331471">~PETScMatrix</a> (void)</td></tr>
<tr class="memdesc:a838b1147e06cd9f9eb3795437d331471"><td class="mdescLeft">&#160;</td><td class="mdescRight">Destructor. <br /></td></tr>
<tr class="separator:a838b1147e06cd9f9eb3795437d331471"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a59a581e26a983fbd21d2accce938bc08"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a59a581e26a983fbd21d2accce938bc08"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a59a581e26a983fbd21d2accce938bc08">Identity</a> ()</td></tr>
<tr class="memdesc:a59a581e26a983fbd21d2accce938bc08"><td class="mdescLeft">&#160;</td><td class="mdescRight">Define matrix as identity matrix. <br /></td></tr>
<tr class="separator:a59a581e26a983fbd21d2accce938bc08"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a1391f8c21e85ef8a89765a7979ea404f"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a1391f8c21e85ef8a89765a7979ea404f"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a1391f8c21e85ef8a89765a7979ea404f">Diagonal</a> ()</td></tr>
<tr class="memdesc:a1391f8c21e85ef8a89765a7979ea404f"><td class="mdescLeft">&#160;</td><td class="mdescRight">Define matrix as a diagonal one. <br /></td></tr>
<tr class="separator:a1391f8c21e85ef8a89765a7979ea404f"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ab2dc6d9a8fed22c57529176cd4643804"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ab2dc6d9a8fed22c57529176cd4643804"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ab2dc6d9a8fed22c57529176cd4643804">Diagonal</a> (const T_ &amp;a)</td></tr>
<tr class="memdesc:ab2dc6d9a8fed22c57529176cd4643804"><td class="mdescLeft">&#160;</td><td class="mdescRight">Define matrix as a diagonal one with diagonal entries equal to <code>a</code> <br /></td></tr>
<tr class="separator:ab2dc6d9a8fed22c57529176cd4643804"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a16ceb7ad1b3c1fafb09a7cc8bc0895c0"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a16ceb7ad1b3c1fafb09a7cc8bc0895c0">setAIJ</a> (const vector&lt; int &gt; &amp;nnz)</td></tr>
<tr class="separator:a16ceb7ad1b3c1fafb09a7cc8bc0895c0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a87177f92fc184e0b042e6810bc509384"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a87177f92fc184e0b042e6810bc509384">setAIJ_MPI</a> (const vector&lt; int &gt; &amp;diag_nnz, const vector&lt; int &gt; &amp;off_nnz)</td></tr>
<tr class="separator:a87177f92fc184e0b042e6810bc509384"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a05b7680ed419a3422be7b7989f307496"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a05b7680ed419a3422be7b7989f307496">setMesh</a> (<a class="el" href="classOFELI_1_1Mesh.html">Mesh</a> &amp;mesh, size_t dof=0)</td></tr>
<tr class="memdesc:a05b7680ed419a3422be7b7989f307496"><td class="mdescLeft">&#160;</td><td class="mdescRight">Determine mesh graph and initialize matrix.  <a href="#a05b7680ed419a3422be7b7989f307496">More...</a><br /></td></tr>
<tr class="separator:a05b7680ed419a3422be7b7989f307496"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ad55999723924ca14c1364e1924b26980"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ad55999723924ca14c1364e1924b26980">setPartition</a> (<a class="el" href="classOFELI_1_1Partition.html">Partition</a> &amp;p)</td></tr>
<tr class="memdesc:ad55999723924ca14c1364e1924b26980"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set a <a class="el" href="classOFELI_1_1Partition.html" title="To partition a finite element mesh into balanced submeshes. ">Partition</a> instance in the class.  <a href="#ad55999723924ca14c1364e1924b26980">More...</a><br /></td></tr>
<tr class="separator:ad55999723924ca14c1364e1924b26980"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:af5dd9b12edc2df7325fe8d3d948383d2"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#af5dd9b12edc2df7325fe8d3d948383d2">setRank</a> (int np, int r=0)</td></tr>
<tr class="memdesc:af5dd9b12edc2df7325fe8d3d948383d2"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set number of processors and processor rank.  <a href="#af5dd9b12edc2df7325fe8d3d948383d2">More...</a><br /></td></tr>
<tr class="separator:af5dd9b12edc2df7325fe8d3d948383d2"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ac1d3aaea63e250685e774e6fb65ab5a0"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ac1d3aaea63e250685e774e6fb65ab5a0"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ac1d3aaea63e250685e774e6fb65ab5a0">setOneDOF</a> ()</td></tr>
<tr class="memdesc:ac1d3aaea63e250685e774e6fb65ab5a0"><td class="mdescLeft">&#160;</td><td class="mdescRight">Activate 1-DOF per node option. <br /></td></tr>
<tr class="separator:ac1d3aaea63e250685e774e6fb65ab5a0"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:af9449887e1b8ca6b5184f015285eacf2"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="af9449887e1b8ca6b5184f015285eacf2"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#af9449887e1b8ca6b5184f015285eacf2">setSides</a> ()</td></tr>
<tr class="memdesc:af9449887e1b8ca6b5184f015285eacf2"><td class="mdescLeft">&#160;</td><td class="mdescRight">Activate Sides option. <br /></td></tr>
<tr class="separator:af9449887e1b8ca6b5184f015285eacf2"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ad739c2f0b2ae98856d555e7c54915c4a"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ad739c2f0b2ae98856d555e7c54915c4a"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ad739c2f0b2ae98856d555e7c54915c4a">setSymmetric</a> ()</td></tr>
<tr class="memdesc:ad739c2f0b2ae98856d555e7c54915c4a"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set matrix as a symmetric one. <br /></td></tr>
<tr class="separator:ad739c2f0b2ae98856d555e7c54915c4a"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a4b7a6c037a0d510ed1977c20bf758bcb"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a4b7a6c037a0d510ed1977c20bf758bcb">DiagPrescribe</a> (<a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;<a class="el" href="group__OFELI.html#gaf5cd2eaaceda19772fd904539819158c">b</a>, const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;u)</td></tr>
<tr class="memdesc:a4b7a6c037a0d510ed1977c20bf758bcb"><td class="mdescLeft">&#160;</td><td class="mdescRight">Impose by a diagonal method an essential boundary condition using the <a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance provided by the constructor.  <a href="#a4b7a6c037a0d510ed1977c20bf758bcb">More...</a><br /></td></tr>
<tr class="separator:a4b7a6c037a0d510ed1977c20bf758bcb"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a97deb13018c1b55b18bfe25a9e339e8d"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a97deb13018c1b55b18bfe25a9e339e8d">setSize</a> (size_t size)</td></tr>
<tr class="memdesc:a97deb13018c1b55b18bfe25a9e339e8d"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set size of matrix (case where it's a square matrix).  <a href="#a97deb13018c1b55b18bfe25a9e339e8d">More...</a><br /></td></tr>
<tr class="separator:a97deb13018c1b55b18bfe25a9e339e8d"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a6cd98f8e60c54091bc2249e99677f8c3"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a6cd98f8e60c54091bc2249e99677f8c3">setSize</a> (size_t nr, size_t nc)</td></tr>
<tr class="memdesc:a6cd98f8e60c54091bc2249e99677f8c3"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set size (number of rows) of matrix.  <a href="#a6cd98f8e60c54091bc2249e99677f8c3">More...</a><br /></td></tr>
<tr class="separator:a6cd98f8e60c54091bc2249e99677f8c3"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a82ca541dbe3735a35549e951e63bec12"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a82ca541dbe3735a35549e951e63bec12">getRange</a> (int istart, int iend)</td></tr>
<tr class="memdesc:a82ca541dbe3735a35549e951e63bec12"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return the range of matrix rows owned by this processor.  <a href="#a82ca541dbe3735a35549e951e63bec12">More...</a><br /></td></tr>
<tr class="separator:a82ca541dbe3735a35549e951e63bec12"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a3ec0cb8e215315aef5111f925f869f8a"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a3ec0cb8e215315aef5111f925f869f8a">setGraph</a> (const vector&lt; std::pair&lt; size_t, size_t &gt; &gt; &amp;I, int opt=1)</td></tr>
<tr class="memdesc:a3ec0cb8e215315aef5111f925f869f8a"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set graph of matrix by giving a vector of its nonzero entries.  <a href="#a3ec0cb8e215315aef5111f925f869f8a">More...</a><br /></td></tr>
<tr class="separator:a3ec0cb8e215315aef5111f925f869f8a"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:af4d11677b77c61b1c2530c0295d9a966"><td class="memItemLeft" align="right" valign="top">T_&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#af4d11677b77c61b1c2530c0295d9a966">operator()</a> (size_t i, size_t j) const </td></tr>
<tr class="memdesc:af4d11677b77c61b1c2530c0295d9a966"><td class="mdescLeft">&#160;</td><td class="mdescRight">Operator <code>()</code>  <a href="#af4d11677b77c61b1c2530c0295d9a966">More...</a><br /></td></tr>
<tr class="separator:af4d11677b77c61b1c2530c0295d9a966"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ac64de7025d9f5712288b3093ca5408b1"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ac64de7025d9f5712288b3093ca5408b1"></a>
size_t&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ac64de7025d9f5712288b3093ca5408b1">getNbRows</a> () const </td></tr>
<tr class="memdesc:ac64de7025d9f5712288b3093ca5408b1"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return number of matrix rows. <br /></td></tr>
<tr class="separator:ac64de7025d9f5712288b3093ca5408b1"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a879356293a4c2120a7f87bf231ce341c"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a879356293a4c2120a7f87bf231ce341c"></a>
size_t&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a879356293a4c2120a7f87bf231ce341c">getNbColumns</a> () const </td></tr>
<tr class="memdesc:a879356293a4c2120a7f87bf231ce341c"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return number of matrix columns. <br /></td></tr>
<tr class="separator:a879356293a4c2120a7f87bf231ce341c"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a68d4c55c333f8db2cc206b90b1e53041"><td class="memItemLeft" align="right" valign="top">size_t&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a68d4c55c333f8db2cc206b90b1e53041">getLength</a> () const </td></tr>
<tr class="memdesc:a68d4c55c333f8db2cc206b90b1e53041"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return length of matrix.  <a href="#a68d4c55c333f8db2cc206b90b1e53041">More...</a><br /></td></tr>
<tr class="separator:a68d4c55c333f8db2cc206b90b1e53041"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ab50a62ca3636fe05fc9b8a2d346bd1b8"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ab50a62ca3636fe05fc9b8a2d346bd1b8"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ab50a62ca3636fe05fc9b8a2d346bd1b8">getMesh</a> (<a class="el" href="classOFELI_1_1Mesh.html">Mesh</a> &amp;mesh)</td></tr>
<tr class="memdesc:ab50a62ca3636fe05fc9b8a2d346bd1b8"><td class="mdescLeft">&#160;</td><td class="mdescRight">Get <a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance whose reference will be stored in current instance of <a class="el" href="classOFELI_1_1PETScMatrix.html" title="To handle matrices in sparse storage format using the Petsc library. ">PETScMatrix</a>. <br /></td></tr>
<tr class="separator:ab50a62ca3636fe05fc9b8a2d346bd1b8"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:af08ec58f7731ea03dd3c5315a72dcdd3"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#af08ec58f7731ea03dd3c5315a72dcdd3">Mult</a> (const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;x, <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;y) const </td></tr>
<tr class="memdesc:af08ec58f7731ea03dd3c5315a72dcdd3"><td class="mdescLeft">&#160;</td><td class="mdescRight">Multiply matrix by vector and save in another one.  <a href="#af08ec58f7731ea03dd3c5315a72dcdd3">More...</a><br /></td></tr>
<tr class="separator:af08ec58f7731ea03dd3c5315a72dcdd3"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a68c86273fa4b307d704b8c9b2e8c6697"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a68c86273fa4b307d704b8c9b2e8c6697">MultAdd</a> (const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;x, <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;y) const </td></tr>
<tr class="memdesc:a68c86273fa4b307d704b8c9b2e8c6697"><td class="mdescLeft">&#160;</td><td class="mdescRight">Multiply matrix by vector <code>x</code> and add to <code>y</code>.  <a href="#a68c86273fa4b307d704b8c9b2e8c6697">More...</a><br /></td></tr>
<tr class="separator:a68c86273fa4b307d704b8c9b2e8c6697"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a9b02aa416b20d307aef3ef2ba164ecb4"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a9b02aa416b20d307aef3ef2ba164ecb4">MultAdd</a> (T_ a, const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;x, <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;y) const </td></tr>
<tr class="memdesc:a9b02aa416b20d307aef3ef2ba164ecb4"><td class="mdescLeft">&#160;</td><td class="mdescRight">Multiply matrix by vector <code>a*x</code> and add to <code>y</code>.  <a href="#a9b02aa416b20d307aef3ef2ba164ecb4">More...</a><br /></td></tr>
<tr class="separator:a9b02aa416b20d307aef3ef2ba164ecb4"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ab561d0b22c12bad47f6d7b95237593ab"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ab561d0b22c12bad47f6d7b95237593ab">set</a> (size_t i, size_t j, const T_ &amp;a)</td></tr>
<tr class="memdesc:ab561d0b22c12bad47f6d7b95237593ab"><td class="mdescLeft">&#160;</td><td class="mdescRight">Assign a value to an entry of the matrix.  <a href="#ab561d0b22c12bad47f6d7b95237593ab">More...</a><br /></td></tr>
<tr class="separator:ab561d0b22c12bad47f6d7b95237593ab"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ab5b687846a29319bfb8577eec442438c"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ab5b687846a29319bfb8577eec442438c">add</a> (size_t i, size_t j, const T_ &amp;a)</td></tr>
<tr class="memdesc:ab5b687846a29319bfb8577eec442438c"><td class="mdescLeft">&#160;</td><td class="mdescRight">Add a value to an entry of the matrix.  <a href="#ab5b687846a29319bfb8577eec442438c">More...</a><br /></td></tr>
<tr class="separator:ab5b687846a29319bfb8577eec442438c"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a5f99e1b114e6f63b520a95628d81c176"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a5f99e1b114e6f63b520a95628d81c176">set</a> (vector&lt; int &gt; &amp;ir, vector&lt; int &gt; &amp;ic, vector&lt; T_ &gt; &amp;val)</td></tr>
<tr class="memdesc:a5f99e1b114e6f63b520a95628d81c176"><td class="mdescLeft">&#160;</td><td class="mdescRight">Assign values to a portion of the matrix.  <a href="#a5f99e1b114e6f63b520a95628d81c176">More...</a><br /></td></tr>
<tr class="separator:a5f99e1b114e6f63b520a95628d81c176"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a8dc08d8cbd80a17fbea793fa2ca1be7f"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a8dc08d8cbd80a17fbea793fa2ca1be7f">operator=</a> (const T_ &amp;a)</td></tr>
<tr class="memdesc:a8dc08d8cbd80a17fbea793fa2ca1be7f"><td class="mdescLeft">&#160;</td><td class="mdescRight">Operator <code>=</code>  <a href="#a8dc08d8cbd80a17fbea793fa2ca1be7f">More...</a><br /></td></tr>
<tr class="separator:a8dc08d8cbd80a17fbea793fa2ca1be7f"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ac8bb3912a3ce86b15842e79d0b421204"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ac8bb3912a3ce86b15842e79d0b421204"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ac8bb3912a3ce86b15842e79d0b421204">clear</a> ()</td></tr>
<tr class="memdesc:ac8bb3912a3ce86b15842e79d0b421204"><td class="mdescLeft">&#160;</td><td class="mdescRight">Set all matrix entries to zero. <br /></td></tr>
<tr class="separator:ac8bb3912a3ce86b15842e79d0b421204"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a5b3c9703591fef88517e4cad780ae64f"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a5b3c9703591fef88517e4cad780ae64f">Laplace1D</a> (<a class="el" href="group__Util.html#gaccfeb6b1e8cf41731fde610549bee67c">real_t</a> h, bool mpi=false)</td></tr>
<tr class="memdesc:a5b3c9703591fef88517e4cad780ae64f"><td class="mdescLeft">&#160;</td><td class="mdescRight">Sets the matrix as the one for the Laplace equation in 1-D.  <a href="#a5b3c9703591fef88517e4cad780ae64f">More...</a><br /></td></tr>
<tr class="separator:a5b3c9703591fef88517e4cad780ae64f"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ada04f559c6aef86159de41938510ca3f"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ada04f559c6aef86159de41938510ca3f">Laplace2D</a> (size_t nx, size_t ny, bool mpi=false)</td></tr>
<tr class="memdesc:ada04f559c6aef86159de41938510ca3f"><td class="mdescLeft">&#160;</td><td class="mdescRight">Sets the matrix as the one for the Laplace equation in 2-D.  <a href="#ada04f559c6aef86159de41938510ca3f">More...</a><br /></td></tr>
<tr class="separator:ada04f559c6aef86159de41938510ca3f"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a462297334738e1512d94504b98b80820"><td class="memItemLeft" align="right" valign="top">int&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a462297334738e1512d94504b98b80820">solve</a> (<a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;<a class="el" href="group__OFELI.html#gaf5cd2eaaceda19772fd904539819158c">b</a>)</td></tr>
<tr class="memdesc:a462297334738e1512d94504b98b80820"><td class="mdescLeft">&#160;</td><td class="mdescRight">Solve the linear system of equations.  <a href="#a462297334738e1512d94504b98b80820">More...</a><br /></td></tr>
<tr class="separator:a462297334738e1512d94504b98b80820"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ad2c1b8da41be364c1351d683f9086bc5"><td class="memItemLeft" align="right" valign="top">int&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ad2c1b8da41be364c1351d683f9086bc5">solve</a> (const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;<a class="el" href="group__OFELI.html#gaf5cd2eaaceda19772fd904539819158c">b</a>, <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;x)</td></tr>
<tr class="memdesc:ad2c1b8da41be364c1351d683f9086bc5"><td class="mdescLeft">&#160;</td><td class="mdescRight">Solve the linear system of equations.  <a href="#ad2c1b8da41be364c1351d683f9086bc5">More...</a><br /></td></tr>
<tr class="separator:ad2c1b8da41be364c1351d683f9086bc5"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a2213133b882c0d42b2b99517eb06b5e4"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a2213133b882c0d42b2b99517eb06b5e4">setSolver</a> (string solver, string prec, <a class="el" href="group__Util.html#gaccfeb6b1e8cf41731fde610549bee67c">real_t</a> toler=1.e-12, int max_it=1000)</td></tr>
<tr class="memdesc:a2213133b882c0d42b2b99517eb06b5e4"><td class="mdescLeft">&#160;</td><td class="mdescRight">Choose solver and preconditioner for an iterative procedure.  <a href="#a2213133b882c0d42b2b99517eb06b5e4">More...</a><br /></td></tr>
<tr class="separator:a2213133b882c0d42b2b99517eb06b5e4"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ae06e66a8f26e28bdf1e0209f7d0d96bb"><td class="memItemLeft" align="right" valign="top">T_ *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ae06e66a8f26e28bdf1e0209f7d0d96bb">get</a> () const </td></tr>
<tr class="memdesc:ae06e66a8f26e28bdf1e0209f7d0d96bb"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return C-Array.  <a href="#ae06e66a8f26e28bdf1e0209f7d0d96bb">More...</a><br /></td></tr>
<tr class="separator:ae06e66a8f26e28bdf1e0209f7d0d96bb"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a7a7b1d2724593a6ec1212247d7964b8e"><td class="memItemLeft" align="right" valign="top">T_&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a7a7b1d2724593a6ec1212247d7964b8e">get</a> (size_t i, size_t j) const </td></tr>
<tr class="memdesc:a7a7b1d2724593a6ec1212247d7964b8e"><td class="mdescLeft">&#160;</td><td class="mdescRight">Return entry <code>(i,j)</code> of matrix if this one is stored, <code>0</code> otherwise.  <a href="#a7a7b1d2724593a6ec1212247d7964b8e">More...</a><br /></td></tr>
<tr class="separator:a7a7b1d2724593a6ec1212247d7964b8e"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a8a4c449305862cc2f96866ebf8f3e1a6"><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a8a4c449305862cc2f96866ebf8f3e1a6">operator Mat</a> () const </td></tr>
<tr class="memdesc:a8a4c449305862cc2f96866ebf8f3e1a6"><td class="mdescLeft">&#160;</td><td class="mdescRight">Casting operator.  <a href="#a8a4c449305862cc2f96866ebf8f3e1a6">More...</a><br /></td></tr>
<tr class="separator:a8a4c449305862cc2f96866ebf8f3e1a6"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a9947b8ce043b9ec45b33fe93c92e2572"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a9947b8ce043b9ec45b33fe93c92e2572"></a>
PetscReal&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a9947b8ce043b9ec45b33fe93c92e2572">getNorm1</a> () const </td></tr>
<tr class="memdesc:a9947b8ce043b9ec45b33fe93c92e2572"><td class="mdescLeft">&#160;</td><td class="mdescRight">Get 1-norm of matrix. <br /></td></tr>
<tr class="separator:a9947b8ce043b9ec45b33fe93c92e2572"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a7a148febf4384f3b2384fae5c9a467b9"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a7a148febf4384f3b2384fae5c9a467b9"></a>
PetscReal&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a7a148febf4384f3b2384fae5c9a467b9">getFrobeniusNorm</a> () const </td></tr>
<tr class="memdesc:a7a148febf4384f3b2384fae5c9a467b9"><td class="mdescLeft">&#160;</td><td class="mdescRight">Get Frobenius norm of matrix. <br /></td></tr>
<tr class="separator:a7a148febf4384f3b2384fae5c9a467b9"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:ab54387d779d3d7204050e51291329f14"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="ab54387d779d3d7204050e51291329f14"></a>
PetscReal&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#ab54387d779d3d7204050e51291329f14">getNormMax</a> () const </td></tr>
<tr class="memdesc:ab54387d779d3d7204050e51291329f14"><td class="mdescLeft">&#160;</td><td class="mdescRight">Get infinity norm of matrix. <br /></td></tr>
<tr class="separator:ab54387d779d3d7204050e51291329f14"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:aae8f063011ae37198a6c96ad00da6188"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#aae8f063011ae37198a6c96ad00da6188">Assembly</a> (const <a class="el" href="classOFELI_1_1Element.html">Element</a> &amp;el, T_ *a)</td></tr>
<tr class="memdesc:aae8f063011ae37198a6c96ad00da6188"><td class="mdescLeft">&#160;</td><td class="mdescRight">Assembly of element matrix into global matrix.  <a href="#aae8f063011ae37198a6c96ad00da6188">More...</a><br /></td></tr>
<tr class="separator:aae8f063011ae37198a6c96ad00da6188"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a377c96d9edc5455d8b81d7f2e75a7437"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a377c96d9edc5455d8b81d7f2e75a7437">Assembly</a> (const <a class="el" href="classOFELI_1_1Side.html">Side</a> &amp;sd, T_ *a)</td></tr>
<tr class="memdesc:a377c96d9edc5455d8b81d7f2e75a7437"><td class="mdescLeft">&#160;</td><td class="mdescRight">Assembly of side matrix into global matrix.  <a href="#a377c96d9edc5455d8b81d7f2e75a7437">More...</a><br /></td></tr>
<tr class="separator:a377c96d9edc5455d8b81d7f2e75a7437"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a34b17b28c9bf2d38a9a0cc02058ccaaa"><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a34b17b28c9bf2d38a9a0cc02058ccaaa">setAssembly</a> ()</td></tr>
<tr class="memdesc:a34b17b28c9bf2d38a9a0cc02058ccaaa"><td class="mdescLeft">&#160;</td><td class="mdescRight"><a class="el" href="classOFELI_1_1Matrix.html" title="Virtual class to handle matrices for all storage formats. ">Matrix</a> assembly.  <a href="#a34b17b28c9bf2d38a9a0cc02058ccaaa">More...</a><br /></td></tr>
<tr class="separator:a34b17b28c9bf2d38a9a0cc02058ccaaa"><td class="memSeparator" colspan="2">&#160;</td></tr>
<tr class="memitem:a8834efe36631472e72fde5b34ce68b35"><td class="memItemLeft" align="right" valign="top"><a class="anchor" id="a8834efe36631472e72fde5b34ce68b35"></a>
void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="classOFELI_1_1PETScMatrix.html#a8834efe36631472e72fde5b34ce68b35">setMPI</a> ()</td></tr>
<tr class="memdesc:a8834efe36631472e72fde5b34ce68b35"><td class="mdescLeft">&#160;</td><td class="mdescRight">Activate MPI option. <br /></td></tr>
<tr class="separator:a8834efe36631472e72fde5b34ce68b35"><td class="memSeparator" colspan="2">&#160;</td></tr>
</table>
<a name="details" id="details"></a><h2 class="groupheader">Detailed Description</h2>
<div class="textblock"><h3>template&lt;class T_&gt;<br />
class OFELI::PETScMatrix&lt; T_ &gt;</h3>

<p>To handle matrices in sparse storage format using the Petsc library. </p>
<dl class="section warning"><dt>Warning</dt><dd>This class is available only when <a class="el" href="namespaceOFELI.html" title="A namespace to group all library classes, functions, ... ">OFELI</a> has been installed with Petsc.</dd></dl>
<dl class="tparams"><dt>Template Parameters</dt><dd>
  <table class="tparams">
    <tr><td class="paramname">T_</td><td>Data type (double, float, complex&lt;double&gt;, ...) </td></tr>
  </table>
  </dd>
</dl>
</div><h2 class="groupheader">Constructor &amp; Destructor Documentation</h2>
<a class="anchor" id="a9708182d1255298f5eb71c68d7270625"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Default constructor. </p>
<p>Initialize a zero-dimension matrix </p>

</div>
</div>
<a class="anchor" id="ac4143fbd170eb4ee9ddb38c0c0b53614"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>nr</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>nc</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Constructor that initializes current instance as a dense matrix. </p>
<p>Normally, for a dense matrix this is not the right class. </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">nr</td><td>Number of matrix rows. </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">nc</td><td>Number of matrix columns. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="af521d98af518239d3362748df7d02ab6"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>size</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Constructor that initializes current instance as a dense matrix. </p>
<p>Normally, for a dense matrix this is not the right class. </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">size</td><td>Number of matrix rows (and columns). </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a91e7a1b127844bd9b2846db0c928d4a9"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1Mesh.html">Mesh</a> &amp;&#160;</td>
          <td class="paramname"><em>mesh</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>dof</em> = <code>0</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Constructor using a <a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">mesh</td><td><a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance from which matrix graph is extracted. </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">dof</td><td>Option parameter, with default value <code>0</code>.<br />
<code>dof=1</code> means that only one degree of freedom for each node (or element or side) is taken to determine matrix structure. The value <code>dof=0</code> means that matrix structure is determined using all DOFs. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a5ec47afc10afd3a677382bbbe4f06819"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="classOFELI_1_1PETScMatrix.html">PETScMatrix</a> </td>
          <td>(</td>
          <td class="paramtype">const vector&lt; std::pair&lt; size_t, size_t &gt; &gt; &amp;&#160;</td>
          <td class="paramname"><em>I</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>opt</em> = <code>1</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Constructor for a square matrix using non zero row and column indices. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">I</td><td>Vector containing pairs of row and column indices </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">opt</td><td>Flag indicating if vectors I is cleaned and ordered (opt=1) or not (opt=0). In the latter case, this vector can have the same contents more than once and are not necessarily ordered </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<h2 class="groupheader">Member Function Documentation</h2>
<a class="anchor" id="a16ceb7ad1b3c1fafb09a7cc8bc0895c0"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setAIJ </td>
          <td>(</td>
          <td class="paramtype">const vector&lt; int &gt; &amp;&#160;</td>
          <td class="paramname"><em>nnz</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">nnz</td><td></td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a87177f92fc184e0b042e6810bc509384"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setAIJ_MPI </td>
          <td>(</td>
          <td class="paramtype">const vector&lt; int &gt; &amp;&#160;</td>
          <td class="paramname"><em>diag_nnz</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const vector&lt; int &gt; &amp;&#160;</td>
          <td class="paramname"><em>off_nnz</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">diag_nnz</td><td></td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">off_nnz</td><td></td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a05b7680ed419a3422be7b7989f307496"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setMesh </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1Mesh.html">Mesh</a> &amp;&#160;</td>
          <td class="paramname"><em>mesh</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>dof</em> = <code>0</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Determine mesh graph and initialize matrix. </p>
<p>This member function is called by constructor with the same arguments </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">mesh</td><td><a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance for which matrix graph is determined. </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">dof</td><td>Option parameter, with default value <code>0</code>.<br />
<code>dof=1</code> means that only one degree of freedom for each node (or element or side) is taken to determine matrix structure. The value <code>dof=0</code> means that matrix structure is determined using all DOFs. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="ad55999723924ca14c1364e1924b26980"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setPartition </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1Partition.html">Partition</a> &amp;&#160;</td>
          <td class="paramname"><em>p</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Set a <a class="el" href="classOFELI_1_1Partition.html" title="To partition a finite element mesh into balanced submeshes. ">Partition</a> instance in the class. </p>
<p>This member function is to be used when parallel computing is considered. </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">p</td><td>Reference to <a class="el" href="classOFELI_1_1Partition.html" title="To partition a finite element mesh into balanced submeshes. ">Partition</a> instance </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="af5dd9b12edc2df7325fe8d3d948383d2"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setRank </td>
          <td>(</td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>np</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>r</em> = <code>0</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Set number of processors and processor rank. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">np</td><td>Total number of processors. </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">r</td><td>Rank of current processor [Default: <code>0</code> </td></tr>
  </table>
  </dd>
</dl>
<dl class="section warning"><dt>Warning</dt><dd>If this member function is not called, only one processor is used and then sequential computing is involved. </dd></dl>

</div>
</div>
<a class="anchor" id="a4b7a6c037a0d510ed1977c20bf758bcb"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void DiagPrescribe </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>b</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>u</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Impose by a diagonal method an essential boundary condition using the <a class="el" href="classOFELI_1_1Mesh.html" title="To store and manipulate finite element meshes. ">Mesh</a> instance provided by the constructor. </p>
<p>This member function modifies diagonal terms in matrix and terms in vector that correspond to degrees of freedom with nonzero code in order to impose a boundary condition. The penalty parameter is defined by default equal to 1.e20. It can be modified by member function <b>setPenal</b>(..). </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in,out]</td><td class="paramname">b</td><td><a class="el" href="classOFELI_1_1PETScVect.html" title="To handle general purpose vectors using Petsc. ">PETScVect</a> instance that contains right-hand side. </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">u</td><td><a class="el" href="classOFELI_1_1PETScVect.html" title="To handle general purpose vectors using Petsc. ">PETScVect</a> instance that conatins imposed valued at DOFs where they are to be imposed. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a97deb13018c1b55b18bfe25a9e339e8d"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setSize </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>size</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Set size of matrix (case where it's a square matrix). </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">size</td><td>Number of rows and columns. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a6cd98f8e60c54091bc2249e99677f8c3"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setSize </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>nr</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>nc</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Set size (number of rows) of matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">nr</td><td>Number of rows </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">nc</td><td>Number of columns </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a82ca541dbe3735a35549e951e63bec12"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void getRange </td>
          <td>(</td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>istart</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>iend</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Return the range of matrix rows owned by this processor. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[out]</td><td class="paramname">istart</td><td>Index of the first local row </td></tr>
    <tr><td class="paramdir">[out]</td><td class="paramname">iend</td><td>Index of the last local row </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a3ec0cb8e215315aef5111f925f869f8a"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setGraph </td>
          <td>(</td>
          <td class="paramtype">const vector&lt; std::pair&lt; size_t, size_t &gt; &gt; &amp;&#160;</td>
          <td class="paramname"><em>I</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>opt</em> = <code>1</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Set graph of matrix by giving a vector of its nonzero entries. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">I</td><td>Vector containing pairs of row and column indices </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">opt</td><td>Flag indicating if vector <code>I</code> is cleaned and ordered (<code>opt=1</code>: default) or not (<code>opt=0</code>). In the latter case, this vector can have the same contents more than once and are not necessarily ordered </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="af4d11677b77c61b1c2530c0295d9a966"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">T_ operator() </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>j</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Operator <code>()</code> </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">i</td><td>Row index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">j</td><td>Column index </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a68d4c55c333f8db2cc206b90b1e53041"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">size_t getLength </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Return length of matrix. </p>
<p>The length is the total number of stored elements in the matrix </p>

</div>
</div>
<a class="anchor" id="af08ec58f7731ea03dd3c5315a72dcdd3"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void Mult </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>x</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>y</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Multiply matrix by vector and save in another one. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">x</td><td>Vector to multiply by matrix </td></tr>
    <tr><td class="paramdir">[out]</td><td class="paramname">y</td><td>Vector that contains on output the result. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a68c86273fa4b307d704b8c9b2e8c6697"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void MultAdd </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>x</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>y</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Multiply matrix by vector <code>x</code> and add to <code>y</code>. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">x</td><td>Vector to multiply by matrix </td></tr>
    <tr><td class="paramdir">[out]</td><td class="paramname">y</td><td>Vector to add to the result. <code>y</code> contains on output the result. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a9b02aa416b20d307aef3ef2ba164ecb4"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void MultAdd </td>
          <td>(</td>
          <td class="paramtype">T_&#160;</td>
          <td class="paramname"><em>a</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>x</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>y</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Multiply matrix by vector <code>a*x</code> and add to <code>y</code>. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">a</td><td>Constant to multiply by matrix </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">x</td><td>Vector to multiply by matrix </td></tr>
    <tr><td class="paramdir">[out]</td><td class="paramname">y</td><td>Vector to add to the result. <code>y</code> contains on output the result. </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="ab561d0b22c12bad47f6d7b95237593ab"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void set </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>j</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const T_ &amp;&#160;</td>
          <td class="paramname"><em>a</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Assign a value to an entry of the matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">i</td><td>Row index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">j</td><td>Column index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">a</td><td>Value to assign to <code>a(i,j)</code> </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="ab5b687846a29319bfb8577eec442438c"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void add </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>j</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const T_ &amp;&#160;</td>
          <td class="paramname"><em>a</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Add a value to an entry of the matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">i</td><td>Row index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">j</td><td>Column index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">a</td><td>Constant value to add to <code>a(i,j)</code> </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a5f99e1b114e6f63b520a95628d81c176"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void set </td>
          <td>(</td>
          <td class="paramtype">vector&lt; int &gt; &amp;&#160;</td>
          <td class="paramname"><em>ir</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">vector&lt; int &gt; &amp;&#160;</td>
          <td class="paramname"><em>ic</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">vector&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>val</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Assign values to a portion of the matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">ir</td><td>Vector of row indexes to assign (instance of class <code>vector</code>) </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">ic</td><td>Vector of column indexes to assign (instance of class <code>vector</code>) </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">val</td><td>Vector of values to assign (instance of class <code>vector</code>) </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a8dc08d8cbd80a17fbea793fa2ca1be7f"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void operator= </td>
          <td>(</td>
          <td class="paramtype">const T_ &amp;&#160;</td>
          <td class="paramname"><em>a</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Operator <code>=</code> </p>
<p>Assign constant value <code>a</code> to matrix diagonal entries </p>

</div>
</div>
<a class="anchor" id="a5b3c9703591fef88517e4cad780ae64f"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void Laplace1D </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="group__Util.html#gaccfeb6b1e8cf41731fde610549bee67c">real_t</a>&#160;</td>
          <td class="paramname"><em>h</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">bool&#160;</td>
          <td class="paramname"><em>mpi</em> = <code>false</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Sets the matrix as the one for the Laplace equation in 1-D. </p>
<p>The matrix is initialized as the one resulting from P<sub>1</sub> finite element discretization of the classical elliptic operator -u'' = f with homogeneous Dirichlet boundary conditions </p><dl class="section remark"><dt>Remarks</dt><dd>This function is available for real valued matrices only. </dd></dl>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">h</td><td>Mesh size (assumed constant) </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">mpi</td><td>true if MPI is used for parallel computing, false if not (sequential), [Default: <code>false</code>] </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="ada04f559c6aef86159de41938510ca3f"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void Laplace2D </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>nx</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>ny</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">bool&#160;</td>
          <td class="paramname"><em>mpi</em> = <code>false</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Sets the matrix as the one for the Laplace equation in 2-D. </p>
<p>The matrix is initialized as the one resulting from P<sub>1</sub> finite element discretization of the classical elliptic operator -Delta u = f with homogeneous Dirichlet boundary conditions </p><dl class="section remark"><dt>Remarks</dt><dd>This function is available for real valued matrices only. </dd></dl>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">nx</td><td>Number of unknowns in the <code>x</code>-direction </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">ny</td><td>Number of unknowns in the <code>y</code>-direction </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">mpi</td><td>true if MPI is used for parallel computing, false if not (sequential), [Default: <code>false</code>] </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a462297334738e1512d94504b98b80820"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">int solve </td>
          <td>(</td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>b</em></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Solve the linear system of equations. </p>
<p>The default parameters are: the Conjugate Gradient method and the Jacobi method for preconditioner. To change these values, call function setSolver before this function </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in,out]</td><td class="paramname">b</td><td>Vector that contains right-hand side on input and solution on output</td></tr>
  </table>
  </dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>Number of actual performed iterations </dd></dl>

</div>
</div>
<a class="anchor" id="ad2c1b8da41be364c1351d683f9086bc5"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">int solve </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>b</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype"><a class="el" href="classOFELI_1_1PETScVect.html">PETScVect</a>&lt; T_ &gt; &amp;&#160;</td>
          <td class="paramname"><em>x</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Solve the linear system of equations. </p>
<p>The default parameters are: the Conjugate Gradient method and the Jacobi method for preconditioner. To change these values, call function setSolver before this function </p><dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">b</td><td>Vector that contains right-hand side </td></tr>
    <tr><td class="paramdir">[out]</td><td class="paramname">x</td><td>Vector that contains the obtained solution </td></tr>
  </table>
  </dd>
</dl>
<dl class="section return"><dt>Returns</dt><dd>Number of actual performed iterations </dd></dl>

</div>
</div>
<a class="anchor" id="a2213133b882c0d42b2b99517eb06b5e4"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setSolver </td>
          <td>(</td>
          <td class="paramtype">string&#160;</td>
          <td class="paramname"><em>solver</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">string&#160;</td>
          <td class="paramname"><em>prec</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype"><a class="el" href="group__Util.html#gaccfeb6b1e8cf41731fde610549bee67c">real_t</a>&#160;</td>
          <td class="paramname"><em>toler</em> = <code>1.e-12</code>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">int&#160;</td>
          <td class="paramname"><em>max_it</em> = <code>1000</code>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Choose solver and preconditioner for an iterative procedure. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">solver</td><td>Option to choose iterative solver among the macros (see PETSc documentation for more details): <ul>
<li>
<code>KSPRICHARDSON</code>: The Richardson iterative method (Default damping parameter is <code>1.0</code>) </li>
<li>
<code>KSPCHEBYSHEV</code>: The Chebyshev iterative method </li>
<li>
<code>KSPCG</code>: The conjugate gradient method [Default] </li>
<li>
<code>KSPCGNE</code>: The CG method for normal equations (without explicitly forming the product <code>A^TA</code> </li>
<li>
<code>KSPGMRES</code>: The GMRES iterative method (see A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems. Y. Saad and M. H. Schultz, SIAM J. Sci. Stat. Comput. Vo|. 7, No. 3, July 1986, pp. 856-869) </li>
<li>
<code>KSPFGMRES</code>: The Flexible GMRES method (with restart) </li>
<li>
<code>KSPLGMRES</code>: The 'augmented' standard GMRES method where the subspace uses approximations to the error from previous restart cycles </li>
<li>
<code>KSPTCQMR</code>: A variant of QMR (quasi minimal residual) developed by Tony Chan </li>
<li>
<code>KSPBCGS</code>: The BiCGStab (Stabilized version of BiConjugate Gradient Squared) method </li>
<li>
<code>KSPIBCGS</code>: The IBiCGStab (Improved Stabilized version of BiConjugate Gradient Squared) method in an alternative form to have only a single global reduction operation instead of the usual 3 (or 4) </li>
<li>
<code>KSPFBCGS</code>: The flexible BiCGStab method. </li>
<li>
<code>KSPCGS</code>: The CGS (Conjugate Gradient Squared) method </li>
<li>
<code>KSPTFQMR</code>: A transpose free QMR (quasi minimal residual) </li>
<li>
<code>KSPCR</code>: The conjugate residuals method </li>
<li>
<code>KSPLSQR</code>: The LSQR method </li>
<li>
<code>KSPBICG</code>: The Biconjugate gradient method (similar to running the conjugate gradient on the normal equations) </li>
<li>
<code>KSPMINRES</code>: The MINRES (Minimum Residual) method </li>
<li>
<code>KSPSYMMLQ</code>: The SYMMLQ method </li>
<li>
<code>KSPGCR</code>: The Generalized Conjugate Residual method </li>
</ul>
</td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">prec</td><td>Option to choose preconditioner in an enumerated variable <ul>
<li>
<code>PCJACOBI</code>: [Default] Jacobi (<em>i.e.</em> diagonal scaling) preconditioning </li>
<li>
<code>PCBJACOBI</code>: Block Jacobi preconditioning, each block is (approximately) solved with its own KSP object </li>
<li>
<code>PCSOR</code>: (S)SOR (successive over relaxation, Gauss-Seidel) preconditioning </li>
<li>
<code>PCEISENSTAT</code>: An implementation of SSOR (symmetric successive over relaxation, symmetric Gauss-Seidel) preconditioning that incorporates Eisenstat's trick to reduce the amount of computation needed </li>
<li>
<code>PCICC</code>: Incomplete Cholesky factorization preconditioners </li>
<li>
<code>PCILU</code>: Incomplete factorization preconditioners </li>
<li>
<code>PCASM</code>: Use the (restricted) additive Schwarz method, each block is (approximately) solved with its own KSP object </li>
<li>
<code>PCLU</code>: Uses a direct solver, based on LU factorization, as a preconditioner </li>
<li>
<code>PCCHOLESKY</code>: Uses a direct solver, based on Cholesky factorization, as a preconditioner </li>
</ul>
</td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">toler</td><td>Tolerance for convergence [Default: <code>1.e-12</code>] </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">max_it</td><td>Maximum number of allowed iterations [Default: <code>1000</code>] </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="ae06e66a8f26e28bdf1e0209f7d0d96bb"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">T_* get </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Return C-Array. </p>
<p>Non zero terms of matrix is stored row by row. </p>

</div>
</div>
<a class="anchor" id="a7a7b1d2724593a6ec1212247d7964b8e"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">T_ get </td>
          <td>(</td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">size_t&#160;</td>
          <td class="paramname"><em>j</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Return entry <code>(i,j)</code> of matrix if this one is stored, <code>0</code> otherwise. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">i</td><td>Row index </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">j</td><td>Column index </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a8a4c449305862cc2f96866ebf8f3e1a6"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">operator Mat </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const</td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Casting operator. </p>
<p>This member functions enables casting an instance of class <a class="el" href="classOFELI_1_1PETScMatrix.html" title="To handle matrices in sparse storage format using the Petsc library. ">PETScMatrix</a> into the <code>Petsc</code> matrix type <code>Mat</code>. This is useful when one wants to usr any Petsc function that is not available in the wrapper (class <a class="el" href="classOFELI_1_1PETScWrapper.html" title="This class is a wrapper to be used when the library Petsc is installed and used with OFELI...">PETScWrapper</a>) or <a class="el" href="classOFELI_1_1PETScMatrix.html" title="To handle matrices in sparse storage format using the Petsc library. ">PETScMatrix</a>. </p>

</div>
</div>
<a class="anchor" id="aae8f063011ae37198a6c96ad00da6188"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void Assembly </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1Element.html">Element</a> &amp;&#160;</td>
          <td class="paramname"><em>el</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">T_ *&#160;</td>
          <td class="paramname"><em>a</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Assembly of element matrix into global matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">el</td><td>Reference to element instance </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">a</td><td><a class="el" href="classOFELI_1_1Element.html" title="To store and treat finite element geometric information. ">Element</a> matrix as a C-array </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a377c96d9edc5455d8b81d7f2e75a7437"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void Assembly </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="classOFELI_1_1Side.html">Side</a> &amp;&#160;</td>
          <td class="paramname"><em>sd</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">T_ *&#160;</td>
          <td class="paramname"><em>a</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p>Assembly of side matrix into global matrix. </p>
<dl class="params"><dt>Parameters</dt><dd>
  <table class="params">
    <tr><td class="paramdir">[in]</td><td class="paramname">sd</td><td>Reference to side instance </td></tr>
    <tr><td class="paramdir">[in]</td><td class="paramname">a</td><td><a class="el" href="classOFELI_1_1Side.html" title="To store and treat finite element sides (edges in 2-D or faces in 3-D) ">Side</a> matrix as a C-array </td></tr>
  </table>
  </dd>
</dl>

</div>
</div>
<a class="anchor" id="a34b17b28c9bf2d38a9a0cc02058ccaaa"></a>
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void setAssembly </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td></td>
        </tr>
      </table>
</div><div class="memdoc">

<p><a class="el" href="classOFELI_1_1Matrix.html" title="Virtual class to handle matrices for all storage formats. ">Matrix</a> assembly. </p>
<p>This function assembles matrix (begins and ends) </p>

</div>
</div>
</div><!-- contents -->
</div><!-- doc-content -->

</div> <!-- class text -->
</div> <!-- main menu -->
<!-- Google Search -->
<br><br>
<center>
  <div id="cse" style="width: 80%;">Loading</div>
  <script src="http://www.google.com/jsapi" type="text/javascript"></script>
  <script type="text/javascript"> 
     google.load('search', '1', {language : 'en', style : google.loader.themes.SHINY});
     google.setOnLoadCallback(function() {
        var customSearchOptions = {};  var customSearchControl = new google.search.CustomSearchControl(
        '012396140824982761142:-qrywxcfi_o', customSearchOptions);
        customSearchControl.setResultSetSize(google.search.Search.FILTERED_CSE_RESULTSET);
        customSearchControl.draw('cse');
     }, true);
   </script>
   <style type="text/css">
     .gsc-control-cse {
     font-family: Verdana, sans-serif;
     border-color: #DAE0E5;
     background-color: #DAE0E5;
   }
   .gsc-control-cse .gsc-table-result { font-family: Verdana, sans-serif; }
   input.gsc-input { border-color: #B6BEC5; }
   input.gsc-search-button {
     border-color: #B6BEC5;
     background-color: #D0D1D4;
   }
   .gsc-tabHeader.gsc-tabhInactive {
     border-color: #999999;
     background-color: #EEEEEE;
   }
   .gsc-tabHeader.gsc-tabhActive {
     border-color: #999999;
     background-color: #999999;
   }
   .gsc-tabsArea { border-color: #999999; }
   .gsc-webResult.gsc-result,
   .gsc-results .gsc-imageResult {
     border-color: #FFFFFF;
     background-color: #FFFFFF;
   }
   .gsc-webResult.gsc-result:hover,
   .gsc-imageResult:hover {
     border-color: #D2D6DC;
     background-color: #EDEDED;
   }
   .gsc-webResult.gsc-result.gsc-promotion:hover {
     border-color: #D2D6DC;
     background-color: #EDEDED;
   }
   .gs-webResult.gs-result a.gs-title:link,
   .gs-webResult.gs-result a.gs-title:link b,
   .gs-imageResult a.gs-title:link,
   .gs-imageResult a.gs-title:link b { color: #0568CD; }
   .gs-webResult.gs-result a.gs-title:visited,
   .gs-webResult.gs-result a.gs-title:visited b,
   .gs-imageResult a.gs-title:visited,
   .gs-imageResult a.gs-title:visited b { color: #0568CD; }
   .gs-webResult.gs-result a.gs-title:hover,
   .gs-webResult.gs-result a.gs-title:hover b,
   .gs-imageResult a.gs-title:hover,
   .gs-imageResult a.gs-title:hover b { color: #0568CD; }
   .gs-webResult.gs-result a.gs-title:active,
   .gs-webResult.gs-result a.gs-title:active b,
   .gs-imageResult a.gs-title:active,
   .gs-imageResult a.gs-title:active b { color: #0568CD; }
   .gsc-cursor-page { color: #0568CD; }
   a.gsc-trailing-more-results:link { color: #0568CD; }
   .gs-webResult .gs-snippet,
   .gs-imageResult .gs-snippet,
   .gs-fileFormatType { color: #5F6A73; }
   .gs-webResult div.gs-visibleUrl,
   .gs-imageResult div.gs-visibleUrl { color: #5F6A73; }
   .gs-webResult div.gs-visibleUrl-short { color: #5F6A73; }
   .gs-webResult div.gs-visibleUrl-short { display: none; }
   .gs-webResult div.gs-visibleUrl-long { display: block; }
   .gs-promotion div.gs-visibleUrl-short { display: none; }
   .gs-promotion div.gs-visibleUrl-long { display: block; }
   .gsc-cursor-box { border-color: #FFFFFF; }
   .gsc-results .gsc-cursor-box .gsc-cursor-page {
      border-color: #999999;
      background-color: #FFFFFF;
      color: #0568CD;
   }
   .gsc-results .gsc-cursor-box .gsc-cursor-current-page {
      border-color: #999999;
      background-color: #999999;
      color: #0568CD;
   }
   .gsc-webResult.gsc-result.gsc-promotion {
      border-color: #D2D6DC;
      background-color: #D0D1D4;
   }
   .gsc-completion-title { color: #0568CD; }
   .gsc-completion-snippet { color: #5F6A73; }
   .gs-promotion a.gs-title:link,
   .gs-promotion a.gs-title:link *,
   .gs-promotion .gs-snippet a:link { color: #0066CC; }
   .gs-promotion a.gs-title:visited,
   .gs-promotion a.gs-title:visited *,
   .gs-promotion .gs-snippet a:visited { color: #0066CC; }
   .gs-promotion a.gs-title:hover,
   .gs-promotion a.gs-title:hover *,
   .gs-promotion .gs-snippet a:hover { color: #0066CC; }
   .gs-promotion a.gs-title:active,
   .gs-promotion a.gs-title:active *,
   .gs-promotion .gs-snippet a:active { color: #0066CC; }
   .gs-promotion .gs-snippet,
   .gs-promotion .gs-title .gs-promotion-title-right,
   .gs-promotion .gs-title .gs-promotion-title-right *  { color: #333333; }
   .gs-promotion .gs-visibleUrl,
   .gs-promotion .gs-visibleUrl-short { color: #5F6A73; }
   </style>
</center>
<!-- End of Google Search Box -->
<br><br><br></p>
<div id="foot_bar">
  Copyright &copy; 1998-2019 Rachid Touzani&nbsp;
</div>
</div> <!-- main -->
</body>
</HTML>
